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We present a family of spin ladder models which admit exact solution for the ground state 
and exhibit non-Haldane spin liquid properties as predicted recently by Nersesyan and Tsvelik 
[Phys. Rev. Lett. 78, 3939 (1997)], and study their excitation spectrum using a simple variational 
ansatz. The elementary excitation is neither a magnon nor a spinon, but a pair of propagating 
triplet or singlet solitons connecting two spontaneously dimerized ground states. Second-order 
phase transitions separate this phase from the Haldane phase and the rung-dimer phase. 
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It is well known that one-dimensional (Id) Heisenberg 
antiferromagnets can exhibit several types of disordered 
"quantum spin liquid" phases. The spin-i chain has a 
unique disordered gapless ground state with power-law 
decay of spin correlations, and its elementary excitations 
are pairs of spinous carrying spin i ||l[] . The ground state 
of the frustrated S = ^ chain with sufficiently strong 
next-nearest neighbor interaction is doubly degenerate, 
the excitations are also spinon pairs, but the spectrum 
is gapful &-01; in presence of any finite exchange alter- 
nation along the chain the spinon pairs get confined into 
well defined magnon excitations Q. The spin-1 (Hal- 
dane) chain has a unique spin-liquid ground state with a 
gap above it formed by a triplet of magnons carrying spin 
5=11^. The two-leg S — ^ ladder, i.e., two Heisen- 
berg S — ■^ chains coupled by a transverse exchange, 
also has a disordered gapful ground state with magnons 
as elementary excitations , and is believed to be essen- 
tially in the same phase as the Haldane chain, as well as 
frustrated S = ^ chain with alternating exchange [pHlO[ . 

Recently, Nersesyan and Tsvelik mU have proposed an 
interesting example of a Id "non-Haldane spin-liquid" 
which has a gapped spectrum but whose excitations are 
neither spinous nor magnons. Using field-theoretical ar- 
guments, they have shown that under certain conditions 
a two-leg S = ^ Heisenberg ladder with additional leg- 
leg biquadratic interaction enters spontaneously dimer- 
ized phase with the excitation spectrum determined by 
the two-particle continuum, and identified the elemen- 
tary excitations as pairs of singlet and triplet domain 
walls connecting the two dimerized ground states. 

In this Letter we present a set of models which ex- 
hibit non-Haldane spin liquid properties as predicted by 
Nersesyan and Tsvelik, and whose ground state can be 
found exactly. We study their excitation spectrum within 
a simple variational approach, and discuss phase transi- 
tions into the Haldane and other phases. 



We start from a more general ladder Hamiltonian 
which includes also transverse interaction along the lad- 
der diagonals and two additional biquadratic interac- 
tions. The model is described by the Hamiltonian 



H — 2_^ J{^lA^l,i+l + 82,^82, i+l) + Jr^l,iS2,i 



(1) 



+ ^(Sl,iSi_i+l)(S2,iS2,i+l) + Jd{^l.i^2,i+1 + ^2.i^l.i+l) 
+^{(Sl,iS2,i+l)(S2,iSi^i+l) — (Si,iS2,i)(Si_i+iS2,i+l)} , 

where the indices 1 and 2 distinguish lower and upper 
legs, and i labels rungs. The model considered by Ners- 
esyan and Tsvelik corresponds to Jd = K = 0. To 
construct the ground state ^0 for the Hamiltonian (|^), 
we will use the technique of matrix-product (MP) states 
||l2| , |l3| . We start from the following ansatz: 

*o = tr{.gi(ii) • g2iu) ■ ■■g2N-i{u) ■ g2N{u)} , (2) 
g^{u) = u ■ l\s), + (T+^\t+i), + (J-^\t^i)i + cr°|fo). 

Here js), and ji^), are the singlet and triplet states of 
the i-th rung, 2A'' is the total number of rungs (peri- 
odic boundary conditions are assumed), 1 is the 2x2 
unit matrix, a'^ are the Pauli matrices, and u, u are free 
parameters. A simpler version of this ansatz with u = u 
describes several known examples of valence bond type 
(VBS) states, e.g., at w = the wave function ^0 is the 
ground state of the effective Affleck-Kennedy-Lieb-Tasaki 
(AKLT) chain |l5[ whose S = 1 spins are composed from 
pairs of S" = ^ spins of the ladder rungs. 



H — 2_^ S„S„_|_i — /3(S„ • S„^i) 



(3) 



at f3 = — -i, and for u = 1 or u = 00 one obtains two 
degenerate ground states of the Majumdar- Ghosh model 
. Originally (H) with u — u was proposed in Ref . |l4| as 
a variational wave function, and recently it was used by 
us [pO| to construct another class of exact ground states 



for a more general ladder model. In the following, we set 
u ^u, then the state vE'o is dimerized and the translation 
for one rung leads to a different state with the same en- 
ergy. The ansatz (B) obeys rotational symmetry, i.e., ^o 
is a global singlet |1^,16{. 

The Hamiltonian (W can be rewritten as a sum of 
identical local terms coupling neighboring rungs, H = 
^j(/ii^i+i — Eq). Let us demand that the wave function 
(0) is a zero-energy ground state of H (which can always 
be achieved by the appropriate choice of £'o), then the 
following requirements have to be fulfilled [Q : (i) the lo- 
cal Hamiltonian ft.i,i+i has to annihilate ^Qj which, due 
to the product property of (||), means that all elements 
of the two matrix products gi{u)gi+i(u), gi(u)gi+i{u) 
should be zero-energy eigenstates of /ii,i+i; (ii) the other 
eigenstates of /i^.i+i should have positive energy. Those 
two conditions fix the structure of the local Hamiltonian 
as follows: 



J 



hi,i+i ^ 2^ 2^ \j\'ipjM){i^jM\ 

J=0,l,2 Af=-,/ 



(4) 



where the eigenvalues Aj > 0, and \iPjm) are the compo- 
nents of the positive-energy multiplets constructed from 
the states of the four-spin plaquette (i^i + 1): 



iV'oo) = [i + {uu)^]'^^\V3\ss)+uu\tt)j^o}, 
\^,)^[2 + f]-'^^f\tt)j=, + \.st) + \t.s)} 

\i'2) = \tt)j=2, f={u + u)/V2. 



(5) 



Here we use the notation |tt),7=i for the triplet of states 
with the total spin J — 1 constructed, in turn, from two 
triplets on rungs i and i + I, etc. 

Now we demand that the structure (Q) is compati- 
ble with the desired form of the Hamiltonian (yj), which 
yields the connection between the parameters J, Jr, Jd, 
V, K on one hand, and the local eigenvalues \j and 
singlet weight parameters u, u of the ground state wave 
function on the other. Those solutions can be classified 
into the following three types: 

(A) " Checkerboard- dimer" model with K ~ Q, Jd 7^ 0: 

u = ±l, u = Tl, V^U/i, X = 0, (6) 

Ai = 1, Ao = 3a;/8, Aa^ 3(1-2:), < 2: < 1, 
J^ = (8J/3)(2-3a;)/(4-3a;), J^ = ^r/2, J > 0, 



the ground state energy density per rung is Eq 



\J, 



and X is an arbitrary parameter. Two degenerate ground 
states are simply checkerboard-type products of singlet 
bonds along the ladder legs. A generic example from this 
family is the model at a; = | with purely biquadratic in- 
terchain interaction: 
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FIG. Ill The excited states \n)l ^ of the model (0), used in 
Eq. ([L1[). Thick solid lines indicate singlet bonds, and thick 
dashed lines can be either singlets or triplets. Arrows in- 
dicate the "direction" of the singlet bonds [i.e., |si^2) = 

2-1/^(1 Tii2)-UiT2))]. 




FIG. The excitation spectrum of the model (CT). The con- 
tinuum is determined by free two-soliton states, its lowest 
boundary is 16-fold degenerate. The dashed line is deter- 
mined by the Haldane triplet excitation ([13) and indicates a 
variational estimate for bound soliton-antisoliton states. 

This "generic" model lies within the class of Hamiltoni- 
ans considered by Nersesyan and Tsvelik. At x = 1 the 
eigenvalue A2 vanishes, indicating a first-order transition 
into the fully polarized ferromagnetic state. 
(B) Multicritical model: 



u = —u, Ao = 0, A2 = 3Ai, 
Jr^V = 4J/3, Jd = Jr/2, K = 0, 



J>0. 



(8) 



J^ = Ja = K = 0, J= 'SV/4 > . 



(7) 



This model has a remarkable property: any wave function 
^o('") of the form (||) with u = ~u is a ground state with 
the same energy per rung Eq = ~\J ■ One can show that 
two ground state wave functions with different values of 
u are asymptotically orthogonal in thermodynamic limit 
iV ^ 00: (*o(^t)|*o("')) = z^^ z{u,u') < 1, so that the 
degeneracy of the ground state is exponentially large. It 
is easy to observe that the model (@) is a particular case 
of (H) at X = 0, so that the model ^) has another phase 
transition point at x = 0; below we will argue that this 
transition is of the first order. 

(C) Model with two second-order phase boundaries 
with Jd = 0, K y^O: 



u = 


-u, K 


V = 


\o{hu^ 


Ai = 


-- AoISm" 



= J, = Ao(m2-1)(w2+3)/2, j^ = 0, (9) 
+ 2m2 + 9)/4, J = 3Ao(u'* + lOw^ + 5)/16, 
+ 14^2+ l5)/8, A2 == Ao(5m''+ 18^2+ 9)/8, 

the g.s. energy per rung is Eq — — ^Ao(7u'' + 22u^ + 19). 
This is a one-parametric family of models since u is ar- 
bitrary (the parameter Aq just sets the energy scale and 
thus is irrelevant). A particular case u = ±1 again leads 
to the "generic" model (0). One can readily observe 
that at M = or u = cx) the ground state is no more 
dimerized. The state with u = Q describes the ground 
state of an effective S* = 1 chain (||) with /3 = — i; the 
state with m = 00 corresponds to a product of singlet 
bonds on the rungs. It is easy to calculate spin-spin and 
dimer-dimer correlation functions Cs(n) = (5'f^iS'f j^_„) 
and Coin) = (AA+n), here A = Si,i -(81,^+1- Si, i_i): 

Cs{n) - [u^ + 3)-i(z+z_)"(5„,2fc - ^-<5„,2fc+i), (10) 
Coin) = 144uV(w^ + 3)^ z± = {u± lf/{u^ + 3) . 

One can see that the dimer correlations exhibit long- 
range order vanishing for u ^ 0, 00, but remarkably there 
is no exponential tail. The spin correlation length goes 
through zero at u = 1 and diverges at w — > 00, how- 
ever, there is no long-range spin order at u -^ 00 since 
the amplitude of spin correlations vanishes in this limit. 
Thus, one can conclude that the model (C) exhibits two 
second-order phase transitions: into the Haldanc phase 
at u = and into the rung-dimer phase at m = 00. We 
will show below that those transitions are characterized 
by vanishing singlet and triplet gaps, respectively. 

By induction with respect to the ladder length one can 
prove that in cases (A) and (C) the two ground states 
given by the MP ansatz are the only ground states of the 
system. 

Elementary excitations of the model (A) can be easily 
visualized as singlet or triplet diagonal bonds separat- 
ing the two ground states and thus being solitons in the 
dimer order (see Fig. H). Since solitons can be created 
only in pairs, the excitation spectrum is a two-particle 
continuum. To study the scattering soliton states, one 
may consider the ladder with 2N + 1 rungs and periodic 
boundary conditions, and write down a simple singlc- 
soliton variational state with certain value of momentum 
p and parity C = ±1: 



\P) 



C 



E' 



,ip(2n+l)i xC 



(11) 



Here the momenta are defined in terms of the Brillouin 
zone of non-dimcrizcd ladder, so that p G [0, tt]. The 
states \n)jg are shown in Fig. H; in a MP formulation 
they can be written as 

n N 

l")l* " n {92t-l{u) g2t{u)j g'^^nll n 52»(w) .921+1 (u) , 

I— 1 I— n+1 

9l=9{u)-C9{u), gl^^a^g{u) + ^g{u)a^^. (12) 
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FIG. H (a) The gaps of different variational excitations for the 
model (0); (b) the same for model (0). Here Agg, A^j,, and 
A]j denote the gaps of singlet-singlet, triplet-triplet soliton 
pairs, and the Haldane triplet, respectively, and (" = ±1 is 
the parity. 

Here /i = 0, ±1 denotes the z-projection of spin of the 
triplet excitation. Another candidate for the role of the 
elementary excitation is a magnon (the Haldane triplet); 
the corresponding variational wavefunction can be again 
written in the form (\llu with 



N 



\n)H= n [92i~i{u) g2tiu)j g" Yi 92i-i{u) g2t{u) , 

i—1 i—n-\-l 

5Cm ^ Cr^52n-l(u).g2n(u) + Cg2n-l{u)cr^' g2n{u) ■ (13) 

The variational dispersion laws have the following form: 

e{p) =eo/[l + 2coA(z,p)], 

A{z,p) = (cos(2p) - z)/{l + z^ ~2zcos{2p)), (14) 

and for the model (0) one has, in J units, 



4 + 3x 

(4-2C)(4-3x)' 



1/4, 



ZH = 

44- 



^H_ 



0, 



^H_ 



39a; 



6(2-C)(4-3a;) 



(15) 



One can see that the lowest energy is always reached for 
the odd-parity states (C — —1)- The Haldane triplet is 
in this case dispersionless, and has a high energy equal 
to 1. The elementary excitation is a soliton-antisoliton 
pair, and for the scattering states its energy is given by 

E{k, q) = es,t [{k + q)/2] + Es^t [[k - q)/2\ , (16) 

where k and q are the total and relative momentum. For 
a; = I [i.e., for the "generic" model (0) with zero trans- 
verse exchange] the energies of triplet and singlet solitons 
coincide. The lowest boundary E{k) of the continuum 
described by ([ij) at a; = | is shown in Fig. g. The gap 
is given by E{0) — E{tt) — i J, and the lowest excitation 
has a 16-fold degeneracy because the states of a soliton 
pair can be classified into two singlets (ss) and (tt),/=o, 
three triplets (st), (ts), and (ii)j=i and one quintuplet 
{tt)j=2- The energy of the Haldane triplet is lower than 
the continuum boundary in the vicinity of the zone cen- 
ter fc = ^ , indicating possible presence of bound soliton- 
antisoliton states. If the transverse exchange is switched 
on (i.e., a; 7^ |), the singlet-triplet degeneracy is lifted, 
and for x < | (a; > |) the lowest excitation is determined 
by singlet (triplet) solitons, respectively. Behavior of the 
corresponding gaps is shown in Fig. 0(a); one can see 
that for both phase transition points a; = and a; = 1 
the gaps remain finite, which suggests that the transition 
to the "multicritical" state at x = is of the first order. 
The ansatze ([l2|), (|l3| ) can be used for the model (C) as 
well. One again obtains the dispersion laws of the form 
(H), with 

4 = z.(i + c^,-'/')/2, 4 = ^t(i - c^r'^')/2 

eg = 12^V{(«2 + 3)2(1 + C^y')}, (17) 

4 = 4(^2 + 2)/{(,i2 + 3)2(1 -C^y^)}, 



c^ = z^' c4^'(i--«) 



off 



8z. 



1/2 



l/2x 



2(1 + C4') 



.1/2- 



u'-3 



1/2 



(u2 + 3)(i + c4^') 

1/2 _ U^ -I 
^ff ^ „,2 , Q- 



Behavior of the gaps is shown in Fig. ^(b) as a function of 
parameter y — u^/{l + u^). Again, the lowest excitations 
are always soliton pairs. At y ^ the odd-singlet soliton 
gap goes to zero, indicating the second-order transition 
to the Haldane phase. At y — > 1 three gaps (of even- 
singlet and odd-triplet solitons and of the even Haldane 



triplet) vanish, signaling another second-order transition 
into the rung-dimer phase. Actually, it follows from ( [l7| ) 
that at y — > (1) the whole continuum of singlet (triplet) 
soliton pairs collapses to zero. 
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FIG. 1. The states \n)j ^ used in Eq. (h_l|), in a special case of the model (M). Thick solid lines indicate singlet 
bonds, and thick dashed lines can be either singlets or triplets. Arrows indicate the "direction" of the singlet bonds [i.e., 

|S1^2>= 2-1/2(1 Tli2>-UlT2»]. 

FIG. 2. The excitation spectrum of the model (W). The continuum is determined by free two-soliton states, its lowest 
boundary is 16-fold degenerate. The dashed line is determined by the Haldane triplet excitation (hsl) and indicates a variational 
estimate for bound soliton-antisoliton states. 

FIG. 3. (a) The gaps of different variational excitations for the model (|6|); (b) the same for model (|9|). Here Ag^, A^^,, 
and Afj- denote the gaps of singlet-singlet, triplet-triplet soliton pairs, and the Haldane triplet, respectively, and ( = ±1 is the 
parity. 



